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IONIZATION PARAMETER AS A DIAGNOSTIC OF RADIATION AND WIND PRESSURES IN H II REGIONS 

AND STARBURST GALAXIES 
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ABSTRACT 

The ionization parameter U is potentially useful as a tool to measure radiation pressure feedback 
from massive star clusters, as it directly reflects the ratio of radiation to gas pressure and is readily 
derived from mid-infrared line ratios. We consider a number of physical effects which combine to 
determine the apparent value of hi in observations encompassing one or many H II regions. An 
upper limit is set by the compression of gas by radiation pressure, when this is important. The 
pressure of shocked stellar winds and the presence of neutral clumps both tend to reduce hi for a 
given intesity of irradiation. The most intensely irradiated regions are selectively dimmed by internal 
dust absorption of ionizing photons, leading to a bias for observations on galactic scales. We explore 
these effects in analytical and numerical models for dusty H II regions and use them to interpret 
previous observational results. We find that radiation pressure confinement sets the upper limit 
\og w hi ~ — 1 seen in individual regions. Unresolved starbursts are known to display a maximum 
value of ~ —2.3. While lower, this is also consistent with a large portion of their H II regions being 
radiation pressure dominated, given the different technique used to interpret unresolved regions, and 
given the bias caused by dust absorption. We infer that many individual, strongly illuminated regions 
cannot be significantly overpressured by stellar winds, and that even when averaged on galactic scales, 
the shocked wind pressure cannot be large compared to radiation pressure. Therefore, most H II 
regions cannot be adiabatic wind bubbles. Our models imply a metallicity dependence in the physical 
structure and dust attenuation of radiation-dominated regions, both of which should vary strongly 
across a critical metallicity of about one-twentieth solar. 

Subject headings: H II regions — galaxies: starburst — galaxies: general — ISM: bubbles — dust, 
extinction — infrared: galaxies 
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1. INTRODUCTION 

Forbidden line ratios are often employed to determine 
the hardness of the radiation field and the composition 
and physical state of irradiated gas within star-forming 
galaxies. Of the quantities revealed by these ratios, the 

/riH - the ratio of ionizing 
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ionization parameter hi 

photon density to hydrogen density - is of special inter- 
est, as it measures the dimensionless intensity of ionizing 
radiation. Becaus e hi controls the ionization state (e.g., 
iTarter et al.lll969f). it is easily disentangled from other 
parameters (|Dopita et al.ll2000T) ; and because it indicates 
the r atio of ioniz ing radiation pressure to gas pressure 
(e.g.. lKrolildll999l ). it encodes information about the rel- 
ative importance of these forces in star formation feed- 
back. 

It is important, therefore, to understand the phys- 
ical and observational effects which influence our 
estimates of hi in distant galaxies. When de- 
termined on kiloparsec scales, \og w hi appears to 
reach a maximum value of about -2.3 in star- 
burst environments, and indeed ho lds this value 
throughout the inner 500 pc of M 82 dThornlev et all 
20001: iForster Schreiber et all l200ll iCarral et all 11994 
Smith et all 120061: iWestmoauette et all 12007^ These 
points led lThornlev et al.l (j2000l ) to argue that this max- 
imum log 10 hi is a property of the starburst phenomenon. 
However, higher values (up to about -1) are sometimes 
seen in cases where H II regions are individua lly resolved 
(jSniiders et alj|2007t Undebetouw et alJl2009D . 



Why do individual H II regions not exhibit arbitrar- 
ily high ionization parameter, and why are starburst de- 
terminations limited to \og 10 U < —2.3? Both limits 
strongly suggest the operation of a saturation mecha- 
nism. We know of four possibilities: 

Dynamics. The characteristic values of hi could 
originate in the dy namics of expanding H II regions. 
iDopita et al. (2005, hereafter D05), in particular, explain 
the saturation of hi in terms of the 'stalling' of expand- 
ing H II regions as they decelerate to the ambient sound 
speed. Under a specific set of assumptions about the 
dynamics of H II regions and their interaction with am- 
bient gas, this implies a unique value of hi at the time of 
stalling, which imprints itself on the line ratios integrated 
over an entire galaxy. However, there are two problems 
with the assumption s that enter th is expla nation. 

First, D05 follow IWeaver et"al] (|1977D in modeling 
H II regions as though they we re enclosed st ellar wind 
bubbles and therefore use the IWeaver et al.l model to 
describe their motions. Two pieces of evidence indi- 
cate that wind energy escapes H II regions, contrary 
to this assumption. Trapped wind would produce x- 
ray luminosities as much as two orders of magnitude 
higher than what is observed ( e.g., in the Carina Nebula: 
lHarper-Clark fc Murravll2009D . Moreover, D05 find that 
the wind luminosity must be drastically scaled down to 
reproduce the observed range of hi values. (We revisit 
this issue below in £13. H and 

Second, D05 assume that the ambient gas is the warm 
(ionized or neutral) component of the ISM, and there- 
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fore has a unique sound speed of roughly lOkms -1 . Al- 
though this is plausible, many of the ionizing photons 
will be consumed as the dense molecular material or cold 
neutral hydrogen is being cleared away. As the thermal 
sound speed of this matter is quite low, the appropriate 
stalling criterion involves the total ( turbulent, magnetic 
and thermal) a mbient pressure (e.g.. lSlavin fc Coxlll993t 
IMatznedl2002T) . 

Homogeneous mixtures. A separate ex p lanat ion has 
been provided by IForster Schreiber et al] (|2001l . FS01) 
in the context of the central 500 pc starb urst core of 
M82, where their and later observations (|Smith et al.l 
I2006t iWestmoauette et~aTl [2007;) demonstrate that U is 
remarkably uniform. The FS01 model invokes an en- 
semble of neutral clouds with ionized outer layers, scat- 
tered among newly formed stars. The uniformity of IA 
must therefore be sought in the properties of this mix- 
ture, which are not part of the IForster Schreiber et al.1 
model. For this reason we do not investigate mixture 
models further, except to point out that most physical 
mechanisms limiting IA in individual H II regions would 
also operate here. 

Radiation pressure confinement. Another physical 
mechanism fo r the saturation of IA was i ntroduced by 
iBinette et al.l ()1997l ) and elaborated by IDopita et al.1 
(2002, D02) in the context of narrow-line regions in ac- 
tive galactic nuclei. When irradiation is sufficiently in- 
tense, ionized gas develops a pressure gradient to oppose 
the force of photons caught by neutral atoms and dust 
grains, and is squeezed into a thin layer near the ioniza- 
tion front (IF). This prevents gas pressure from becom- 
ing small compared to t he ionizing radiation pr essure 
(as previously argued by iFerland fc Elitzurl [19 84) . The 
corresponding upper limit on log 10 IA depends on the ion- 
izing spectrum and the lines used to infer ionization; a 
rough estimate based on pressure equilibrium gives -1.5. 
While this is significantly above the value of -2.3 seen 
in starbursts, the details of line formation, among other 
effects, can reduce IA to about its observed maximum. 

Dust-limited H II regions. Finally, the observed satu- 
ration could simply be a selection effect, caused by the 
fact that dust absorbs most of the ionizing starlight in re- 
gions with \og 10 U > —2, suppressing their line emission. 
Because extragalactic observations inevitably integrate 
over a wide range of local conditions, the inferred IA will 
be weighted toward regions of lower intensity. Such a 
'dust-limited' state has been shown to successfully re- 
produce the suppression of fine structure lines relative to 
far- infrared luminosity in ultra-luminous infrared gal ax- 
ies (|Voitill992t IBottorff et afl[l998t lAbel et al.ll2009aT) . 

Our goal in this paper is to evaluate these possibili- 
ties and, in the process, to assess whether IA can be a 
useful proxy to measure radiation pressure feedback in 
H II regions in external galaxies. We review the physi- 
cal effects which control the state of ionized gas and the 
values of IA inferred from forbidden line ratios, focusing 
on its saturation in starburst galaxies and the relevance 
of the last two mechanisms. Because radiation pressure 
and dust opacity are both proportional to IA, they occur 
simultaneously. But, since they are truly two separate 
physical mechanisms, we take care not to conflate them. 
In this paper we concentrate on the internal structures 
of H II regions, only evaluating a very simple model for 
their dynamics and galactic populations (<gj). A com- 



panion paper (Verdolini et al. 2012, in prep.) presents a 
new suite of models for the emission from galactic pop- 
ulations of H II regions evolved with more sophisticated 
internal dynamics. 

We begin in § [2] by reviewing the parameters of quasi- 
static, dusty, spherical H II regions. We draw here on 
the models of iDrainel (|2011l hereafter Drll), examining 
in detail their planar, radiation-confined limit (Appendix 
[B]), but we also introduce a new parameter to account for 
an inner, pressurized region of shocked stellar wind. We 
discuss the scales of dust opacity and radiation pressure 
in § 12.11 confinement by radiation pressure and the max- 
imum IA in § 12.21 and the influence of an inner bounding 
pressure in § 12.31 

In §[3] we focus on the extra effects which can reduce the 
maximum value of IA, including stellar winds and neutral 
clumps in the region. Magnetic fields can in principle al- 
low IA to violate its nonmagnetic upper limit, but this 
turns out to be difficult to arrange. We discuss the ef- 
fects of internal and external dust extinctions on limiting 
IA in § [4j In § [5) we provide a survey of IA values obtained 
from modeled MIR line ratios within static and internally 
pressurized dusty H II regions, and demonstrate how stel- 
lar winds act to suppress IA \, S ■ In § El we use our results 
to interpret observations on the scales of individual H II 
regions and entire galaxies. 

2. QUASI-STATIC H II REGIONS: PARAMETERS, 
PROPERTIES, AND SATURATION OF W OB s 

In the static or quasi-static state, H II regio ns have sev- 
eral w ell-known limits: classical H II regions (Stromgrcn 
1939), which lack significa nt radiation pressure or dust 
opaci ty; dusty H II regions (lMathislll971l : lPetrosian etaLl 
Il972h in which radiation pressur e is negligible ; and 
radiation-confined H II regions (|Binette et al.l 119971 : 
IDopita et al.l[20 02) in which both radiation pressure and 
dust opacity are generally important. Recently, Drll has 
thoroughly explored the parameter space of dusty H II 
regions with radiation pressure, demonstrating a three- 
parameter family of models which takes all of these lim- 
its. 

Ionization and force balance within spherical H II re- 
gions are especially simple under the following set of ide- 
alizations (Drll): the ionized gas temperature T, the 
mean energy of ionizing photons = (ft-^)j, the effective 
dust opacity per H atom ad, and the effective recombi- 
nation coefficient a, are all uniform constants; pressure 
from reprocessed radiation, such as the Lyman a line, is 
unimportant; the flow is static, so forces must balance; 
the pressure from stellar winds is neglected (which we 
will relax below in §2.31 and §3.ip ; and there are no addi- 
tional forces from gravity or magnetic fields. 

With these idealizations, H II regions are described by 
seven dimensional parameters: S, the central output of 
ionizing photons per second; Li/c and L n /c, the total 
force of ionizing and non- ionizing photonfl; a, the ef- 
fective recombination coefficient; P/nn = 2.2/csT, the 
thermal pressurcQper H density; PrF, the pressure at the 
ionization front; and o~d- 

1 L n refers to the frequency range where the dust opacity is still 
~ a^; Drll assumes that the two opacities are equal. 

2 Ignoring helium, Drll takes this coefficient to be 2; our value 
of 2.2 assumes He is singly-ionized everywhere. 
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Other dimensional parameters of interest, determined 
by the seven above, include: rt« = P\-p / (2.2k B T), the H 
density at the ionization front; i?st = [3S/ (inanit,)] 1 / 3 , 
the radius of a dust-free Stromgren sphere of uniform 
density n^; n om , the recombination- weighted ion density 
(eq. and 
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127r(2.2/cTc) 2 S"' 



(1) 



the radius of a uniform-density, dust-free Stromgren 
sphere whose gas pressure equals the total, unattenuated 
radiation pressu re L / {Ait R 2 t c) at its boun dary. (Our 
r c h differs from iKrumholz fc Matznerl l2009f s [hereafter 
KM09's] r c h in that it accounts neither for the effects 
of dust extinction, nor for the additional pressure from 
stellar winds.) Yet another is the ionization front 
radius. 

We focus on three dimensionless ratios which control 
the structure and appearance of the H II region^ First 
of these is the characteristic ratio of ionizing radiation 
pressure to gas pressure in a classical Stromgren sphere, 



u 
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0.93 

7^1.53 
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3 , S — lO 51 ^! s 1 , and we have 
cm 3 s _1 ()Storev fc Hummer] 



where ni*=10 4 n^*^ cm 
used a = 2.54 x 10" 13 T, 
1995). Here and elsewhere, we use the subscript 'St 
to denote conditions in a classical, uniform Stromgren 
sphere without attenuation, rather than the actual H II 
region. 

An alternative to T foregoes i?st as a dimensional scale 
in favor of the actual ionization front radius -Rif; 



* = 



Rif 

fch 



(3) 



This is useful if i?ip is observed or derived from the dy- 
namics of the region (KM09); we use "J to parameterize 
our simulation output in § [5] Bear in mind that ^ de- 
creases as radiation pressure becomes more important, 
whereas T increases. 

Our second parameter is the dust opacity in a classical 
H II region, 

Td,st = Cd«-i*-Rst- (4) 

Both r and r^st are proportional to (Sn^) 1 ^ 3 ; as a result 
t c i oc hi in uniform H II regions. 
The third parameter is the dust discriminant 
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2.2k b T <j d c T d: st 



: 16.3(7,; -21^ 



a 3r 
18 13.6eV 



(5) 



21 ~- This sets the relative 



where ad = 10 _ '"CTd,-2i cm", 
importance of dust and neutral H atoms in absorb- 
ing ionizing photons when radiation pressure confines 
the gas: in a situation where gas pressure matches the 

3 A fourth, the total number of H atoms in the classical 
Stromgren sphere, is unimportant. 



ionizing radiation pressure, the dust absorbs a fraction 
7/(1 + 7) of all ionizing photons. Fiducial parameters 
(o-i,_2i = 1,74 = 0.8, = 20 eV) give 7 = 7.4. 

Additional parameters are required to capture the 
spectral form of the incident starlight and dust opac- 
ity, the most relevant of which, /3, compares the radi- 
ation force of non-ionizing starlight to that of ionizing 
starlight; see Drll and Appendix [B] 

For emission in a given line, and therefore for the de- 
termination of line ratios and electron temperatures, the 
ratio of the characteristic density m* to the transition's 
critical density n cr it is a significant additional parameter. 
So long as this ratio is small for all of the lines which are 
either observed or contribute to the collisional cooling of 
the region, the density parameter will not be important. 
We return to this point in <J5] 

Although the above idealization of H II regions is rea- 
sonable and admits a simple parameter space, we wish to 
relax its approximations and assumptions when consid- 
ering observational constraints. Microphysical assump- 
tions, such as the constancy of T, can be relaxed through 
the use of a modern photoionization code like Cloudy or 
MAPPINGS; we take this approach in § [5] Macrophys- 
ical issues, such as stellar winds and bulk flow, must be 
addressed separately (§©. 

2.1. Ionization parameter, radiation pressure, and dust 

opacity 

We are especially interested in observational measures 
of the ionization parameter U(r) = riji/nu- This is di- 
rectly related to the ionizing radiation-to-gas pressure 
ratio 



Si(r) 



-Prad,i(r) 
^gas(r) 



= U(r) 



2.2k B T 



(6) 



These are both functions of position r, and are each 
observable only in an integral sense within the region. 
Whereas U(r) and the ionizing spectrum set the local 
ionization structure and recombination spectrum, any 
recombination line ratio will depend on an average of 
IA weighted by the emissivities of the species involved. 
The value W bs inferred from such line ratios will natu- 
rally reflect this weighting. One observationally relevant 
proxy is therefore the recombination- weighted average 



J n 2 H UdV 



(7) 



where both integrals are over the entire H II region. An- 
other commonly-used proxy is a 'geometrical' estimate 
which relies on an observational determination of the 
characteristic ionized hydrogen density: 
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S 



4irR 2 F n iy0hs c 



(8) 



We also define [E item , s iigcom ]= (ei/2.2k B T)[U em ,U eeom \ 
for consistency. In a classical H II region (uniform den- 
sity, negligible dust absorption, radius i?st), 



= -F. 



(9) 



and, in such a region, Sj 



r by definition. 
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Table 1 

Definitions 



Symbol 



Definition 



Reference 



Note 



P,Bi](r) 
«(r) 

Wgcom 
fobs 
^obs t max 
(^obs)gal 
Wcm 

r 

n 

ftrap,w 
f iv 
/ion 

7 





Ratio of [total, ionizing] radiation pressure to gas pressure Eq. 10 
Ionization parameter n 7 ,i/n 

Geometric estimate S / (Air R^n Q ^ s c) Eq. JSJ 
Observational estimate of W for individual H II region or line of sight 
Maximum value of W b s in a quasi-static, unmagnctizcd layer 

Obs. estimate of U for beams spanning many H II regions Eq. 12311 

Recombination- weighted average of U Eq. J7 1 

Char, radius for radiation-gas pressure equality Eq. Q I 

Radiation pressure parameter, Rrp/^ch Eq. lf3l 

Ionizing radiation pressure parameter, VP" 1 / 2 Li/ L Eq. J2| 

Wind pressure parameter, Pi n V| n /(PiF Vjp — P n Vi n ) Eq. 1131 1 

Wind pressure parameter, 4ir R? n Pi n c / L KM09 
Wind force parameter M w v w c/ L 

Fraction of ionizing photons absorbed by gas rather than dust Drll 

Dust optical depth to ion. photons from source to ionization front Eq. Q 

Dust discriminant, oc cr d T/(eia) Eq. ifSJ, Drll 

Effective ratio of non-ionizing to ionizing radiation force Dr ll 

Recombination-weighted average density Eq. 12611 



Position dependent 
Position dependent 

Tracer dependent 
Tracer dependnet 
Biased by dust 
Theoretical proxy for IA \, S 
Similar to r c b in KM09 
* > 1 implies P rad < P 
T < 1 implies P rad < P 
Q > 1 implies P rad < P 



T(j > 1: dust-limited HII region 

7 > 1 implies > 1 when rad-confincd 

See note in Table [3] 

Proxy for observed density 



These emission-weighted averages are directly related 
to the dust opacity, as the following argument shows. 
Following Drll, we define the local ionizing photon den- 
sity to be S(j>(r)/(4irr 2 c), and the integrated recombina- 
tion rate to be / an 2 dV = f ion S, where / ion = 1 - / dust 
is the fraction of ionizing radiation absorbed by hydro- 
gen rather than dust. Rewriting both the numerator and 
denominator of equation ([7]) , we find 



l 2.2k B T 7 / k 



-Td 



where 



/ <j)(r)<T d n H dr 
J a d n H dr 



(10) 



(11) 



is a dust absorption- weighted average of 4>{r). 

Equation (|10[) is remarkable, because it relates a re- 
gion's dust optical depth to its ionizing photon-to-gas 
pressure ratio in terms of the dust discriminant 7, and 
the order-unity parameters /i on and <f> take definite forms 
in each of the limiting regimes of static H II regions. In 
the context of radiation-confined H II regions, which we 
now discuss, 7 determines the importance of dust ab- 
sorption. 

2.2. Radiation confinement and saturation ofU b s 

In a classical H II region (r <C 1,^ >• 1 ), except 
for a small region near the ionizing source (|Mathewsl 
I1967D , radiation pressure is negligible and has no practi- 
cal effect on the distribution of ionized gas. In this case 
the recombination-weighted radiation pressure parame- 
ter 2 em = (1 + /3)Si !em is much less than unity, as can be 
seen from equation ©. 

This ceases to be true as the illumination becomes 
more intense. All of the momentum in ionizing starlight, 
and (thanks to dust) some of the momentum in non- 
ionizing starlight, is transferred to the ionized gas. In 
a quasi-static region without other forces, the radiation 
force must be balanced by a gas pressure gradient. Many 
have pointed out that when radiation pressure is signif- 
icant (r ^> I,*? <C 1), i onized gas can be effect ively 
excluded from the interior (|Elitzur fc Fer land 198(f) and 



confined to a thin shell near the IF (|Binette et al.lli"997t 
iDopita et all 120021: lDraindl20TTI ). Within this shell the 
total pressure is nearly constant, as we show in Appendix 
[A"l so the rise in gas pressure matches the drop in ra- 
diation pressure. Although the ratio of radiation and 
gas pressures varies throughout the region, its character- 
istic value records a state of pressure balance - either 
S^em — 1, if dust absorption is negligible (7 < 1), or 
S om ~ 1 if it is not (7 > 1). We explore the radiation- 
confined limit more quantitatively in Appendix iBl 

Because S em reaches its maximum value of around 
unity in radiation-confined layers, the recombination- 
weighted ionization parameter U cm does as well (Drll): 
in a dusty layer with S em ~ 1, we expect U cm = 
2.2fc B 7 1 Si iCm /[(l + /?)e 4 ] ~ lCT 1 - 5 . (This estimate rises 
to 1CP 1 in low-dust regions with 7 < 1, because of the 
factor of S/Sj = 1+0.) Moreover, the structure of the 
ionized zone becomes nearly independent of the strength 
of the illumination, apart from an overall scale, so long 
as T > 1. 

Any ratio of forbidden lines, which are all well above 
or well below their critical densities, therefore satu- 
rates iat_a_Jmitevahae in this radiation-confined limit 
(e.g. lBinette et al.1 ll997). The ionization parameter W bs 
which would be inferred from a specific set of line ratios 
tends to be close, but not equal, to U QU1 . Because the ra- 
diation force causes U(r) to be strongly stratified, there 
exists a diffuse inner region where lines preferentially 
form that require a higher ionization state ([Dopita et al.l 
2002). Conversely, low- ionization transitions form in 
outer regions of lower Li. The observationally-inferred 
value U bs therefore limits to a maximum value which 
depends on the lines being observed in addition to all 
the parameters apart from the strength of illumination: 
the chemical abundances, dust properties, and incident 
spectrum. In other words 



(12) 



where W bs. max depends on the lines in question and these 
other parameters, but not on the radiation intensity (un- 
less radiation pressure leads the density to cross one of 
its critical values). Likewise, if density-dependent tracers 
are used to derive ni :0 bs and U„ com , then this will exhibit 
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an upper limit in the vicinity of max(W cm ) ~ 10 -1 ' 5 in 
the radiation-confined state. In §[5] we demonstrate the 
saturation of U a h s for a selection of mid-infrared emission 
lines with high critical densities using the Cloudy code. 
We find in Figure [3] that U b s derived from a particular 
[Ne III] / [Ne II] line ratio exhibits a maximum at 10~ 197 
for our fiducial parameters. 

Because of the effects just described, the saturation of 
emission line ratios is a feature of the radiation-confined 
state in which photon momentum dominates over gas 
pressure gradients in the structure and dynamics of H II 
regions. As we shall see in 21 however, a selection effect 
due to dust attenuation (§S|) also modulates the apparent 
value of U bs in large-scale observations of dusty galaxies 
(7 > 1). Both effects are important in the interpretation 
of galactic line ratios. 

In Appendix [51 and Figures [13] and [14j we present the 
properties of H II regions in the radiation-confined state, 
i.e., the planar limit of Drll's equations in the absence of 
a finite pressure at the inner boundary. This limit, which 
corresponds to V ^S> 1 (and 1 in the definitions of 

^2.3[) . is parameterized only by the dust discriminant 7 
and the spectral parameter /3, which quantifies the non- 
ionizing radiation pressure. Although these solutions are 
only a subset of Drll's three-para meter model f amily , 
they improve upon the treatment bv lDopita et al.l (|2003[ ) 
by accounting for the gas density gradient induced by ra- 
diation pressure. They show that: (1) Sj -cm saturates at 
a value which is near unity and only logarithmically de- 
pendent on 7. (2) Dust absorbs most ionizing radiation, 
so Td ~ 1 — 3 for astrophysically relevant values of 7. As 
Drll discusses, 7 ~ 10 at solar metallicity unless dust 
is effectively removed. Therefore radiation-confined H II 
regions are also marginally 'dust-limited'; see ij4.ll (3) 
For relatively dust-free gas with 7 < 1, the planar solu- 
tion is not compact in radius. It must either be bounded 
inward by a bubble of hot gas, or be matched to a spher- 
ical solution (as in Drll)|3 

In the radiation-confined shell, the column density of 
ionized gas takes a characteristic maximum value: cither 
~ 2.2kBTc/(aei) if dust absorption is negligible, or sev- 
eral times a if it is not. What happens if the mean 
column density of the ambient medium is less than four 
times this shell column, so there is not enough matter to 
achieve it? In this case, it is easy to show that the initial 
Stromgren sphere forms outside f c h, so that VP > 1 and 
r < 1: the radiation-dominated phase never occurs. 

2.3. Pressurized inner boundary 

The ionizing sources of H II regions (typically stars, 
star clusters, or AGN) often produce fast outflows which 
decelerate or are deflected before they reach the IF, 
thereby creating regions of hot and pressurized gas within 
the photo-ionized zone. What new parameter should we 
introduce to capture the effect of such zones on the H II 
regions' emission? The ratio of outer to inner pressures is 
not appropriate, because the two are equal in a classical 
H II region and also in a wind-confined shell. Neither is 

4 Along with the relatively short dust drift time scale (Drll), this 
suggests that dust-free gas may escape inward of the dusty layer 
even when 7 is initially large. We do not account for this effect in 
our equations and models, as we assume a uniform mixture of gas 
and dust. 
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Figure 1. Parameter space of quasi-static, spherical, un- 
magnctized, dusty H II regions surrounding central wind bub- 
bles. Wind-confined, radiation-confined, dust-limited, and clas- 
sical (Stromgren) regions are indicated. Black contours denote 
the radial dust optical depth Td from the source to the ionization 
front; blue contours indicate the radial thickness of the ionized 
zone, AR/Rip. We use thick contours to separate fiducial regimes. 
For this figure we adopt /3 = 1, 7 = 7.4; the radiation-confined 
state is encompassed by the dust-limited state so long as 7 > 1. 

the volume fraction of the interior zone, because photo- 
ionized gas can be confined to a thin shell by radiation 
pressure just as well. Instead we choose 



P V 

J m v li] 



PifV if - P m V m 



(13) 



where PwVw — Pi n V ln is the difference of the product 
of pressure and volume between the IF and the inner 
edge of the H II region (the outer edge of a pressur- 
ized bubble). This combination reflects the strength of 
the interior pressure, both for classical and radiation- 
dominated H II regions. In fact ft directly reflects the 
contribution of the interior bubble to the total kinetic 
ener gy of the region, which, a ccording to the virial theo- 
rem (jMcKee & Zweibel|[l992T) . controls the expansion of 
it and any shell of neutral matter it pushes outward (e.g., 
KM09). 

In Figure [1] we indicate the physical regimes of spheri- 
cal, quasi-static, dusty H II regions surrounding a central 
wind bubble. 

Quasi-static HII regions typically have an inward gra- 
dient of U. For classical H II regions this is simply be- 
cause the ionizing flux is greater at smaller radii, while 
for radiation-confined regions it is due to the gas den- 
sity gradient. In either case, replacing a region's inner 
portion with hot gas has the effect of removing its most 
highly photo-ionized zone. This effect will be greatest for 
those lines more which form preferentially in regions of 
highest U. See ij3.1l for a discussion of the effect on U b s - 

In radiation-dominated (r > 1) and relatively dust 
free (7 < 1) regions, low-density ionized gas fills the in- 
terior ( Appendix [B]). Because this low-density gas is ideal 
for the formation of lines sensitive to regions of high U, 
we expect these line ratios to be especially sensitive to 
interior pressure when T ^> 1 and 7 < 1. 

3. FACTORS AFFECTING Wobs IN INDIVIDUAL H II 
REGIONS 

Several physical effects modify U b s relative to its theo- 
retical maximum value U Q h s max within an individual H II 
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region. In the following subsections we consider stellar 
winds, neutral clumps, dust absorption of non-ionizing 
light, and magnetic fields. Of these, only magnetic fields 
act to raise U b s relative to W bs,max- Although highly- 
magnetized regions can in principle have U a h s > W bs,max) 
we shall see that this is very difficult to arrange in prac- 
tice. 

3.1. Pressurization by shocked stellar winds 



medium ((Weaver et alj|1977l ). 



As discussed above in § 12 .31 any force which clears ion- 
ized gas from the interior of an H II region will have 
the effect of reducing the characteristic value of U b s , es- 
pecially for those lines which form in regions of highest 
U. The most relevant force is the pressure of a shocked 
wind produced by the same source that ionizes the re- 
gion. Most such regions are ionized by hot stars, which 
inject fast line-driven winds. However the significance 
of stellar wind pressure depends on the wind strength, 
e.g. the ratio f w = M w v w c/L between the momentum 
fluxes of winds and starlight, and the trapping of wind 
energy within the region. KM09 combine these into a 
single parameter f% T &p,w = 47rPi n i?? n c/L > f w which en- 
capsulates the contribution from shocked winds relative 
to the direct momentum flux of starlight. 

So long as the combined force of winds and radia- 
tion compresses the ionized layer into a thin shell, the 
enclosed volume is constant and il reduces a compari- 
son of pressures. Since Prp ~ P; n + L/(47rr 2 c) in this 
case (Appendix [Aj ignoring on the last term the factor 
[1 + (5(1 — e~ Td )]/(l + (5) which accounts for incomplete 
absorption), we expect 



ft ~ 4irr 2 cP in /L = /, 



trap,w- 



(14) 



Further, if we consider a typical location where half of 
the incident radiation has been absorbed and its pressure 
applied to the gas, the same approximations lead to the 
estimate 

' (15) 



2/i 



trap, to 



+ 1' 



although, in the case where confinement is accomplished 
by radiation alone (/trap,™ "C l)i the values of S cm plot- 
ted in Figure [Ml are more accurate. 

We refer to section 3.1 of KM09 for a detailed discus- 
sion of the characteristic values of f w and /trap,™ one 
might expect around massive star clusters. For our pur- 
poses, it suffices to consider three limiting cases: (1) 
perfectly trapped, adiabatic bubbles of shocked w ind, 
as envisioned i n the m odels of iWeaver et al.l (Il977fl and 
IKoo fc McKei (|1992f) and employed by iDopita et al.l 
ill HH), which ha ve f tTm ,„ ~ U,/fe > 1: 
(2) leaky wind bubbles (|Harper-Clark fc Murravl 120091 
KM09), for which /trap,™ ~ 1 in the case of a very mas- 
sive cluster of solar metallicity (KM09); and (3) very 
weak wind bubbles with /trap,™ *C 1, such as one might 
expect if the central source consists of low-metallicity or 
lower-mass main sequence stars. 

In the first case (adiabatic wind bubbles) f2 3> 1 and 
Uohs will be limited to values well below W bs,max- In- 
deed, U seom can easily be c o mput ed from the self-similar 
solutions of IKoo fc McKed ([1992D : using L = (1 + (5)Sei 
and evaluating for the case of a steady wind in a uniform 



2.2k B T t R 1F 
3(1 + P)&i f w v w 



(16) 



The factor of two in the first equality arises be- 
cause a pressure-confined ionized layer is nearly uni- 
form and optically thin to dust; therefore U drops lin- 
early in radius from the value U gcom at its inner edge. 
Since v w ~ 10 3 kms _1 from a young cluster which 
fully samples the IMF, equation ([TB]) yields U gcom ~ 

1 _3 ' 8 /j 1 (PiF/10kms _1 ). Equation (pU assumes the 
ionization front is trapped within the swept-up ambient 
matter, which requires 



S^n m > 1.5 



1 4 - fl IF6 



(1 + /3) 3 / 



w u w8 



20 eV 



(17) 



where uho is the ambient H density in cm 3 , = 
10i?iF6 kms , and v w — 10 3 v w s kms -1 . 

In the second case (leaky bubbles), the wind force is of 
order ftrnp,wL/c and is therefore significant when /trap,™ 
exceeds unity. So long as the ionized layer is thin, the 
suppression of U \, s relative to Wobs.max is of order 3 as 
given by equation (TT51) . 

In the last case (very weak winds), W bs will be unaf- 
fected by their presence, although high-W lines from the 
innermost regions may be (especially when 7 < 1). In 
our modeling effort (§ Ej), we consider a wide range of 
possible values for CI. 

Spherically symme t ric, partially-radiative bubble s 
(|Koo fc McKed [19921: IMatzner fc Krasnopolskvl [2008), 
which have /trap,™ lower than the adiabatic case but well 
above unity, are not especially relevant here because they 
require slower winds than is typical of ionizing stars. 
Fast-wind bubbles can also become radiative, however, 
because of mass entra inment from e vapor ating clumps 
(see EU). Note that iMcKee et ail (|1984f ) argued that 
this effect prevents the filling factor of shocked wind, and 
therefore also /trap,™, from exceeding unity. 

As an aside, we note that some ultracompact H II re- 
gions are probably formed in quasi-static layers confined 
by the ram pressure of protostellar winds, which are both 
highly collimated and much stronger than radiation pres- 
sure. Comb ining the outflow force-luminosity correlation 
re ported by I Wu et al.l (|2004[) with the collimation model 
of lMatzner fc McKed (|1999f ). 



/™ =s 7.4 



(10 4 L Q /L) 
sm(8) 2 + 10- 



().:S5 



(18) 



for an angle 9 from the wind axis and a bolometric 
source luminosity L. Such regions are typically in the 
first, wind-dominated regime, even well away from the 
wind axis. This ram-pressure confined phase can only 
exist in directions where the wind itself is fully ionized. 
The alternate cas e, in which the IF for ms in the wind, 
was considered bv lTan fc McKed ([20031 ). It is inertially- 
confincd in the sense of ij3.4l below, and can exhibit 

£4,bs > Z^obs .max ■ 

3.2. Clumps within H II regions 

H II regions are often porous and clumpy instead of 
spherical, as we invoked above in the context of wind 
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trapping. Here we wish to point out that dumpiness 
reduces the inferred U h s relative to a spherical region of 
the same radius. The net effect is therefore to reduce 
Uohs relative to W b s ,max, although the upper limit does 
not necessarily change. 

In the absence of significant radiation pressure, a 
clump subjected to ionizing radiation is first crushed b y 
the recoil of its photo-evaporative flow (|Bertoldil 1989). 
If the clump is dense enough to survive the first implo- 
sion phase, it comes i nto a quasi-static equilib rium with 
the induced pressure (jBertoldi fe McKedl 199(1 . in which 
its radius of curvature is r c (on the symmetry axis, or the 
point most directly illuminated by starlight). Ionization 
balance in the presence of this accelerating and diverging 
photoionized flow implies 



S 



ujan^Tr 



where rii C is the ionized gas density just outside the clump 
IF at this point, R c is the distance from the clump to 
the c entral source, and uj = 0.1 — 0.2 (|Bertoldi fc McKeel 
1990) The parameter which describes the characteristic 
ionizing radiation-to-gas pressure ratio (Hj) for this out- 
flow is therefore 



/ ujaR c S 



2.2kTc V 4?rr 2 



1/2 



(19) 



Compared to conditions at a spherical IF of radius -Rif, 



T 2 C , r c R w 

1 IF n c 



(20) 



A clump generally has a smaller value of T, and therefore 
U, than its parent region, because its radius of curvature 
r c is significantly less than the ionization front radius Rw 
of the region. 

A clump's photoevaporation flow ends at a decelera- 
tion shock caused by the pressure from ionized gas in 
the parent region, stellar winds, or radiation. As T c in- 
creases toward unity, radiation pressure requires that the 
shock approach the sonic point of the photoevaporation 
flow, and when T c > 1, radiation impedes the outflow of 
ionized gas. For r c S> 1, the ionized gas is confined by 
radiation pressure to a thin layer on the surface of the 
clump, leading to a saturation of Hj C m and lA oos , exactly 
as described in § !2.2l and AppendixlBl Just a s before, dust 
opacit y will be significant so long as 7 > 1. IDopita et al.l 
(|2002ft sketch the flow for T c < 1. 

However there is one significant difference between the 
effects of dumpiness and stellar winds. The ionized layer 
of a clump cannot be hydrostatic: it must either evapo- 
rate away, or be dragged along the clump surface by the 
tangential component of the radiation momentum flux, 
and also possibly by the stress from a turbulent boundary 
layer between photoionized clump gas and shocked stellar 
winds. Such wind and photon-driven flows are a plausi- 
ble source of supersonic line widths in giant H II regions, 
for photoevaporative flows only produ ce motions of up to 
abou t twice the ionized sound speed (jBertoldi fc McKeel 
1990). It is not at all clear that a dynamical equilibrium 
exists which would prevent a clump with r c > 1 from 
being dynamically disrupt ed, although a magnetic buffer 
zone may allow this (e.g.. lDursll2007D . 



3.3. Magnetic fields 

Magnetic fields supply an additional form of pressure 
which is available to balance the force of radiation, and 
therefore have the potential to raise U oos above W bs,max- 
However this is somewhat difficult to arrange. 

Consider a planar shock of speed v s running into a 
medium of hydrogen density no, magnetic field per- 
pendicular to the shock normal B±q and Alfven speed 
vao = Sj_o/(47r/ifl-no) 1 / 2 . The post-shocked fluid is in 
contact with (and driven forward by the pressure of) an 
ionized layer. Ignoring any change in pressure across the 
post-shocked layer as well as streaming of the ionized gas 
(both of which can be order- unity effects), the density of 
ionized gas at the contact discontinuity is (v s /vao) 2 tiq if 
it is pressure-supported, and \> r 2(v s /vAo) n o if it is mag- 
netically supported. Magnetic pressure support therefore 
requires 

v s v A0 > v^c 2 (21) 



lOkms 



is the ionized sound 



where Cj 1 
speed. Magnetic support is therefore only 

relevant for shock speeds in excess of about 
35(n /100cm- 3 ) 1 / 2 T 4! 4(30^G/B ± o)kms- 1 . 

Magnetic fields should be relatively more important 
within clumpy regions, because of the tendency for 
streaming and turbulent flows to amplify them toward 
equipartition with the kinetic energy. 

3.4. Inertially- Confined and Impulsive Regions 

Our discussion so far has focused on regions which are 
in force balance because ionized gas resides within them 
for long enough to be crossed by sound waves. However, 
some regions will be too young (or have crossing times 
too long) to develop into this state, such as the regions 
created in cosmological reionization. In others, ionized 
gas flows through the region at speeds well in excess of 
the ionized sound speed; these include ionized winds and 
jets, ionized ejecta, ionized zones of supersonic turbu- 
lence, and H II r egions around ru naway stars. Indeed, 
champagne flows (jWhitw orth 1979|) are density-bounded 
to one side and are only approximately in force balance 
on the side toward the ionization front. None of these 
impulsive or inertially-confined flows is limited to 5 < 1 
and U \y S < ^obs.max for the reasons we have outlined: 
therefore, higher values of U bs could be taken as spec- 
tral evidence of this state. 

Loosely speaking, an observation of U oos above 
W bs,max implies that either a non-thermal form of pres- 
sure, such as ram pressure or magnetic fields, balances 
radiation pressure, that another force (such as gravity) 
is at play, or that there has not been time to establish 
equilibrium. 

3.5. Additional Effects: Recombination Line Pressure 
and Photoelectric Heating 

It is worth noting a couple additional effects involving 
the interaction among light gas, and dust: 

Lyman-a trapping and dust absorption. In the absence 
of dust, photons from the Ly-a line scatter many times 
before escaping. Once the scattered line radiation is suf- 
ficiently strong its pressure will aug ment gas pressure, 
leading to an increase of W b s ,max (See lElitzur fc Ferlandl 



Table 2 

Mid-Infrared Forbidden Lines a 



Species 


A (jim) 


IP (eV) 


n crit (cm 3 ) 


[Aril] 


6.99 


15.76 


2.0 x 10 5 


[Arlll] 


8.99 


27.63 


3.0xl0 5 


[Nell] 


12.81 


21.56 


6.5xl0 5 


[Nelll] 


15.55 


40.96 


1.3x10 s 


[S in]] 


18.71 


23.34 


2.0xl0 4 


[SIV] 


10.51 


34.79 


6.0xl0 4 



1986, for an analysis of line pressure in the AGN con- 
text). The effect is analogous to the role of magnetic 
fields, except that the line pressure is proportional to the 
ionizing radiation pressure. At solar metallicity, a typi- 
cal grain population limits the Ly-a line pressure within 
the ionized zone to 1 only ~ 6% of the ionized gas pressure 
(|Hennev & Arthur! Il998h . but line pressure is significant 
at low metallicity or in regions cleared of dust. Even in 
the absence of dust, the dynamical effect of lin e pressure 
limits logm U 5~, — 1 within qu asi-static regions (|Williamsl 
U97l iFerland fc Elitzurlfl98l . 

Dust photoelectric effect. Photoelectrons ejected from 
dust grains represent an additional sources of heat and 
ionization, while collisions and recombinations with dust 
grains are (less sig nificant) sinks. These processes, con- 
sidere d in detail by Ivan Hoof et al.l (|200l . [Groves et al.l 
(2004), and others, are fully accounted for within the 
Cloudy code. They depend both on the grain abundance 
(captured here by 7) and the grain size distribution. 

Of these, trapped line radiation tends to make W bs,max 
to decrease as 7 increases, whereas photoelectric heating 
has the opposite effect by raising the ionized gas temper- 
ature in dustier gas. 

4. DUST-BOUNDED REGIONS AND SATURATION OF U 
FROM SELECTIVE DUST ATTENUATION 

Emission lines usually arise from regions large enough 
to encompass many local environments, so their ratios re- 
flect an average of the local conditions. Dust absorption 
has the effect of dimming high-W regions, suppressing 
their contribution to the average and causing a saturation 
of the apparent ionization level. This bias arises from the 
dust absorption of ionizing radiation within H II regions, 
and from the absorption of the forbidden lines them- 
selves, predominantly outside the regions themselves. 

4.1. Internal dust extinction 

If a region is sufficiently dusty, then the dust opti- 
cal depth will exceed unity before radiation pressure be- 
comes important. The critical value of the dust discrim- 
inant 7 is apparent if we combine its definition (eq. [5| 
with the characteristic ionizing radiation-to-gas pressure 
ratio in a classical Stromgren sphere (Eq. [9J : 



T"d.st 



47 

T' 



(22) 



Although this refers to the classical case of a filled 
Stromgren sphere, for which both r^st "C 1 and Sj iein <§C 
1, it demonstrates that, when 7 > 3/4, dust extinction 
of ionizing photons is significant even for some regions 
which are not radiation dominated. Conversely, when 



7 < 3/4, radiation pressure confinement can occur when 
dust extinction is not strong. The former scenario (dust- 
limited regions) is consistent with our expectations of 
Milky Way conditions, so long as dust is not cleared 
from the gas; Figure [T] maps the physical regimes in this 
case. The latter applies to low-metallicity environments, 
or more broadly if dust clearing is effective. (We will 
not distinguish the critical 7 from unity in our following 
discussion.) 

Because recombination light is always redder than the 
photons which ionize an H II region, and because dust 
opacity increases with frequency, the line photons of in- 
terest tend to escape the ionized gas relatively unab- 
sorbed. Nevertheless, because dust consumes ionizing 
photons, the rate of recombinations is reduced in dust- 
dominated regions. The ionization efficiency f lon drops 
from unity (for a classical H II region) to an asymp- 
totic value of order (1 + §)Td/ (7 — 1) (in the planar limit, 
when 7 > 1: Appendix [Bj as radiation pressure and dust 
extinction become significant; see Drll's Eq. (21). For 
the parameters of interest, the minimum value of /i on is 
about 0.3. 

Because of this, recombination light from the regions 
of highest U b s is suppressed by the factor / ion - so that, 
on a galactic scale, U Q b s is reduced relative to W bs,max- 
Specifically, so long as the extinction of line photons 
( 34.2[) is negligible, the inferred value of U a b s on galac- 
tic scales will be approximately the emission-weighted 
average of all the regions encompassed by the observa- 
tion: we call this (W bs) ga i- If we label a population of 
Af regions by their ionizing luminosities and apparent 
ionization parameters, then 
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where in the second equality we assume dhf /d]nS oc 
S^^ 1 and define (x) s to be the average of x among re- 
gions whose ionization rates lie in a narrow range around 
S. Galactic H II region luminosity functions obey £ ~ 
at high S, but are significantly shallower at lower lu- 
minosities ( S < 10 50 s~ 1 , where the IMF sampling is 
incomplete; iKennicutt et al.l fl989f) . Each decade of S 
above this break therefore contributes about equally to 
(Wobs)g a i, while those below it are increasingly insignifi- 
cant. 

In general, a galaxy with more massive star clusters (a 
population extending to higher S) and a denser circum- 
cluster medium (so that regions are more compact) will 
have regions characterized by higher values of U bs and 
therefore attain a higher (M bs} sa y However, the decline 
in /ion accompanying the increase in Ta,st suppresses the 
contribution of those regions with highest U Q b s ■ So long 
as the H II region population includes a range of intensi- 
ties and a significant range of and fi on , the net effect 
will be that (I4 bs) sa \ is an underrepresentation of the 
starlight-weighted average of U D b s (the result of taking 
/ion — > 1 in eq. |2"3"]) . In particular, (U Q b s }„ a i will be closer 
to the value it has when r^st = lj 
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It is important to remember that this selection ef- 
fect does not apply to individual H II regions, or parts 
thereof, which saturate at a higher value (W bs,max, ^2.2[) . 
It operates only when starlight photons are distributed 
rather evenly across environments with a range of U Q b s 
and /ion, so that variations of /i on can affect the appar- 
ent value of U bs ■ For this reason, it does not prevent 
(^obs) ga i from limiting £/ bs,max if all the H II regions in 
a galaxy are radiation-confined. 

We demonstrate these points in Figure [2 for which 
we generate a toy model of the H II region distribution 
within a galaxy. We start by examining a p opulation 
of 10 6 associat ions generated for a prior work (jMatznerl 
12001 based on lMcKee fc Williams! I1997L in which £ = 
at high S), extracting from it the luminosity function 
dM/d\nS and the typical ionizing lifetime (ti on ) s as 
functions of S. We then realize a scenario in which all 
H II regions expand into density profiles with p oc r -1 , 
taking the mean column density E = 2rp to be constant 
and the same for all regions. (For simplicity, we neither 
vary E with S, nor alter the upper limit of S with E - 
although both effects are likely to occur in real galaxies.) 
In each luminosity bin we evaluate the time evolution 
^(t) using KM09's equation (13), which accounts for the 
combined force due to radiation, partially trapped winds 
(if any), and photo- ionized gas. We limit regions' growth 
according to their finite initial size and stalling in the 
finite hydrostatic pressure of the environment (KM09's 
§ 2.3). We assume for this model that stellar winds are 
weak: log 10 £1 = —1.5. With this assumption, we can 
construct fi on (S,t) and U b s (S,t) (from the Nelll/Nell 
line ratio) using our own numerical models (S}5|). Taking 
the time average for < t < (ti on ) s to get (/ion) 5 and 
(/ion^obs)s, we arrive at (W bs) Ea i- For comparison we 
also plot the starlight- weighted averages of U b s , /ion, and 
the radiation-to-gas force ratio F ra d/F to t = (1 + vj/ 1 / 2 ) -1 
(see KM09 eq. 10). Although the dust selection effect is 
not strong in our toy model, amounting to a suppression 
of (Wobs)g a i by only ~ 7 — 12% over a relevant range of 
E, we suspect it might be more significant if some of our 
toy model assumptions were relaxed. 

Our toy model does not account for heterogeneity of 
the ambient medium, which permits the existence of a 
diffuse ionized component; it is likely to under-estimate 
(/ion) g a i and over-estimate (Mobs) ga \ for this reason. 

4.2. External dust extinction 

In the expansion of an H II region, most of the ambi- 
ent medium is swept before the ionized gas rather than 
being ionized itself. Therefore, although dust within the 
H II region does not severely attenuate recombination 
line photons, the surrounding matter poses a thicker bar- 
rier. This external dust extinction is also an increasing 
function of U bs, as the following argument shows. 

Suppose that each star cluster, of mass M, forms 
within an overdensity of radial profile p = por~ k , and 
moreover that, statistically speaking, po oc M J . The toy 
model of § 14.11 adopted k — 1 and j — 0, but k could be 
as low as or as high as 2, and j is likely to be positive. 
For the star clusters massive enough to sample the IMF, 
L oc S oc M . 

How then does the dust optical depth of the neutral 
matter outside an H II region scale with U within it? We 



9 



icrf ur] 10? 10' 




E (cga) 

Figure 2. Influence of internal dust absorption on the fraction of 
ionizing photons caught by H atoms rather than dust ({/ion} ga i> top 
panel), galaxy-averaged ratio of radiation to total (gas+radiation) 
force ((F r ad/^tot} ga i, top panel) and apparent ionization parame- 
ter ((^obs)gai, bottom panel), in a toy model for a galactic pop- 
ulation of H II regions ( i)4.1| l. For com parison we also plot the 
starlight-averaged value of W D bs (eq. 1231 with /i on — > 1). In this 
model, the luminosity function and lifetimes of the driving asso- 
ciations are derived from M onte-Carlo modeling by sampling the 
stellar IMF (Matzncr 2002). Regions expan d due to a combination 
of rad iation and photoionized gas pressures (Krumholz & Matzner 
2009) in a medium of definite column density E, and are taken to 
absorb and emit light as in the numerical models of ij5] 

restrict attention to k > 1, so that optical depth varies 
significantly with radius. When the region size is Rif, 
the external dust optical depth at any frequency is then 

Td,ext oc Rifp(Rif) — PoR\F' k i at the same time, the 
radiation pressure parameter oc R\f/S, so that 

r rf ,oxt oc S^ 1 '^ 1 -* (25) 

for massive clusters. When k > 1, regions of a given 
S are most enshrouded when they are most radiation- 
dominated (i.e., when '5 is lowest). The region luminos- 
ity increases or decreases Td at fixed depending on 
whether j > k — 1 or not. However, regions with higher 
S tend to also be more radiation dominated (i.e., have 
lower , J). When j > k — 1 it is clear that Td decreases 
as *f> increases. Insofar as this external dust is optically 
thick at the frequencies of the lines being observed, this 
will tend to suppress emission from high-W b s regions, re- 
ducing this quantity relative to its theoretical maximum 
value. Note that, as we saw in § 12.21 the radiation domi- 
nated state which produces W bs.max is in fact impossible 
if the column density of the ambient medium is too low. 

The external attenuation of line photons is minimized 
by considering the longest possible transition wave- 
lengths; we focus on mid-infrared lines in § \5\ It should 
also be mitigated by holes in the swept-up neutral shell, 
as w ell as dumpiness of the ambient dust distribution 
(e.g., iFischera fc Dopitdl2005l ). 

5. NUMERICAL MODELS OF QUASI-STATIC H II REGIONS 

In this section we present numerical results for spher- 
ical, dusty, quasi-static H II regions. Our goal is to ac- 
count explicitly for the possibility of an inner boundary 
pressure due to shocked wind (as parameterized by Q) 
and to demonstrate the saturation of line ratios and in- 
ferred U hs due to radiation pressure confinement. We 
consider only a few ionizing spectra and do not vary the 
composition of gas or the grain population; a full explo- 
ration of these parameters is beyond the scope of this 
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work. 

We use C l oudy version 08.00, last described by 
iFerland et al.l (|1998l ). to account for many important mi- 
crophysical effects including the photoelectric effect, col- 
lisional cooling, and the pressure due to optically thick 
recombination lines. However we consider only quasi- 
static regions in perfect force balance, and do not ac- 
count for secular effects, such as grain drift, which might 
lead to inhomogeneities in the composition. 

5.1. Spectral synthesis and photoionization models 

Using Starburst99 (|Leitherer et al.lll999h . we gener- 
ate the ionizing continua from coeval star clusters of 
different ages, all of which we assume are massive 
enough to fully sample the stellar initial mass func- 
tion, which we take to have exponents -1.3 and - 
2.3 between stellar mass boundaries 0.1, 0.5, and 120 
Mq. We employ the Geneva high mass-loss evolution- 
ary tracks with solar metallicity. These are optimized 
for modeling at mospheres of high mass s tars and are 
recommended by iMaeder fc Mevnetl (|1994D . We adopt 
Pauldrach/Hillier atmospheres, as these include non- 
LTE and line-blanketing effects (|Smith et al.l 120021 ) for 
O stars dPauldrach et all 120011 ) and Wolf-Rayet stars 
(jHillier fe Milled 1 1998[ )~The combination of the Geneva 
high mass-loss tracks and Pauldrach/Hillier atmosphere 
is recommended when Wol f-Rayet stars are important 
(| Vazquez fc Leithererl [20051. Starburst99 output spec- 
tra are recorded from to 11 Myr with 0.5 Myr steps. 

Starburst99 output continuum spectra are fed into 
Cloudy as the ionizing continuum of each simulated H II 
region. Each H II region is spherical and in perfect 
force balance; we allow radiation pressure to exceed gas 
pressure, in contrast to Cloudy's default setting. We 
adopt Cloudy's default ISM abundances and dust grain 
size distributions. Each calculation stops where tem- 
perature drops to 100 K, and so encapsulates the IF. 
Each set of the simulations outputs the integrated lu- 
minosity of selected MIR emission lines form Table [H 
including [ArIIl] A 9.o^m, [Ar ll]x7.0fim, [NeIIl] A i 5 . 5)UrrL , 
[NeIl] A i2. 8 /zm, [SlIl]Ai8.7/im, and [S Iv] A io.5^m- 

5.2. Mapping line ratios to U a h s 

In both observational and theoretical studies it is nec- 
essary to translate some observed line ratios into a set 
of physical parameters, including U; as we have al- 
ready discussed, we use U b s to designate the apparent 
U which characterizes a single region. A standard defi- 
nition for Uobs is the value of U at the inner boundary of 
an (often uniform) slab which reproduces the observed 
lines; however this has pitfalls both for low and high 
hi (JT}, and ignores the possible role of wind pressur- 
ization. We choose instead to associate U h s with the 
value of U within a homogeneous mixture of gas and 
radiation. To do so, we require a library of the line 
luminosities emitted by such homogeneous mixtures, in 
which physical parameters such as U, n, and input spec- 
trum are varied. We construct this library using the 
very innermost zone of a sequence of Cloudy models. Se- 
lected results are presented in figures [3] to [5] Figure [7] 
demonstrates how well the [Ne III] /[Nell] ratio reflects 
U. Other line ratios, [ArllljAs.gg/zm/IArll^e.gg^m and 

[S IV] A io.5/Ltm/[S M]Ai8.7/im, are a l so most sensitive to the 






4 


4 — 










log 


-3.5 


3.5 


3.5- 




-3 


3 


3— 






— 2.5 — 


2.5 — 





log 10 [Ar III] /[Aril] . 



-2.5— 




0.5 



logio(iff)(cm 3 ) 

Figure 3. Ratio of Arlllg.gg^m to ArIl6.gg Mm luminosities as a 
function of U and nu in homogeneous mixtures of solar metallicity 
gas with the ionizing spectrum of a fully sampled coeval star cluster 
at 2 Myr age. The change in line ratio across the critical densities 
of 2 — 3 X 10 5 cm -3 is clearly visible if not dramatic. 
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Figure 4. As in figure[3] but for Nelllis.ss ^ m /Ne II12.81 nm- 
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Figure 5. As in figure[3] but for S IVio.5i lMm /S III18.71 urn- 
intensity of radiation field, i.e., to U. 

5.3. Model parameters 

We work with a fixed dust population and several in- 
put spectra corresponding to fully sampled, coeval star 
clusters of different ages. Therefore, while the ratio /3 
of non-ionizing to ionizing radiation force varies greatly 
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Table 3 

Coeval Massive Cluster Ionizing Spectra 



Table 4 

Physical parameters of starburst galaxies 



Age (Myr) 


g, (eV) 


Ln/ Li 


/3 a 


loglO °d 


7 





20.6 


0.61 


0.61 


-20.9 


11.2 


2 


19.1 


1.25 


1.12 


-20.8 


13.8 


4 


19.2 


2.81 


2.51 


-20.8 


13.2 


6 


16.7 


16.7 


11.3 


-20.7 


19.3 



a For /3 we list the ratio of non-ionizing to ionizing radiation 
force transferred to the dust grains by an unattenuated in- 
put spectrum, differentiating this from L n /Li : the definition 
used by Drll. 

from young (2 Myr) to relatively old (6 Myr) cluster spec- 
tra, the dust discriminant 7 is almost constant. We list 
the relevant quantities in Table [3] Whenever we are not 
exploring age dependence, we employ spectra from 2 Myr 
clusters, as this is halfway through their typical ionizing 
lifetimes. 

For each spectrum our goal is to map the line luminosi- 
ties, line ratios, and resulting U h s as function of ^ and 
r, i.e., to scan the parameter space displayed in Figure 
[TJ Practically, we accomplish this by varying the cen- 
tral luminosity and innermost density of the region, at 
a fixed inner radius. For this study we restrict ourselves 
to regions with ionized densities well below the critical 
densities for the transitions being considered (Table [2]) 
so that density does not enter as a third parameter. Our 
plots in this section should therefore not be used for re- 
gions with pressures exceeding 10 8 ks cm -3 for [S III] and 
[Siv], or 10 9 ks cm -3 for the others. Moreover we keep 
the metallicity of gas and stars fixed, and vary only the 
age of the stellar population (rather than its IMF or other 
properties). 

In addition to the line luminosities, we compute /i on 
from the ratio of total H/3 emission to the amount of 
H/3 emission expected in the absence of dust. We also 
compute <SVi em , a useful observational diagnostic of the 
strength of ionization, where 



_ Jn 3 dV 



(26) 



the recombination-averaged density - a proxy for the ob- 
servationally inferred density n b s - For regions without 
internal wind bubbles, Sn em determines T and (Sj2j 
Drll); but pressurization breaks this relation. 

5.4. Results 

We begin by comparing the physical structures of a 
radiation-confined shell (log 10 * = —1.09, log 10 f2 = 
— 1.56) with those of a wind-confined shell (log 10 v & = 
1.71, log 10 ft = 0.85) in Figured] As expected, radia- 
tion confinement leads to a dramatic gradient of electron 
density and pressure, whereas wind pressure causes the 
region to be nearly uniform. 

Figures [5] to [TU] show W obs derived from [Ne III] / [Ne II] 
ratios for cluster ages 2, 4, and 6 Myr. Immediately ap- 
parent is the saturation of U b s in the radiation-confined 
state, log 10 ^ < and log 10 H < 0, as is its suppression 
in wind-confined shells (log 10 fi > 0). 

At age 4 Myr, the most massive stars have left the 
main sequence, so the ionizing flux drops and U b s values 
decrease accordingly. The pattern of saturation in the 
radiation-dominated quadrant is, however, unchanged. 



Object 


n e (cm 3 ) 


log 10 5 (s- 1 ) 


^obs 


References 


M82 central 500 pc 
NGC 253 
NGC 3256 
Antennae 


10 - 600 
430±»° 
<1400 


54.09 
53.00 
52.30 - 53.78 


-2.1 - -2.6 
-2.2 - -2.6 

-2.3 

-1.57 


a 

b, c, d 
c 
f 



a IForstc r Schrcibc Fet al.l 1200111 

b ICarral etal| ^19941) 

c lEngelb racht et al] U998I) 

d IThornlev et al.l 1120001) 

e IRov et all t2005l) 

f ISniiders et all 1 120071) 

By 6 Myr, only stars with mass < 30 M Q are still on 
the main sequence, which are responsible for producing 
only 7.5% of the ionizing flux at zero age. Therefore U Q hs 
values decreased gradually comparing to that at younger 
ages. Note that U h s still saturates in the same manner 
as the 2 Myr and 4 Myr cases. 

A second set of contours in the $>-fl plots, labeled with 
dashed black lines, indicate \og 1Q (Sn cm ). This quantity 
increases as \& decreases, which implies that dust opacity 
increases as radiation pressure also becomes more impor- 
tant. This is consistent with the analytical prediction in 
§ [5J Sn cm values are weakly correlated with f2 in the 
log 10 < 1 regimes where stellar winds are negligible. 
Then the contours become a linear function of both * 
and f2 when log 10 f2 > 1. 

As we have discussed in § [4] dust grains absorb a frac- 
tion 1 — /i on of the ionizing photons in H II regions. We 
estimate /; on using the ratio of total H/3 emission to the 
dust-free value H/3o, which we calculate by assuming that 
all ionizing photons are used for photoionization, and 
each ionizing photon is responsible for one recombina- 
tion. The emitted H/3 luminosity - which we tally zone- 
by-zone to avoid extinction of the line photons - reflects 
the suppression of ionizations due to dust absorption of 
starlight. The variation of H/3/H/3o with if? and f2 is 
shown in Figure [TTI with solid blue contours, for the ion- 
izing spectrum of a 2 Myr-old cluster. Again, dashed 
black contours denote log 10 (iS , ?i cm ), providing a calibra- 
tion with respect to previous U b s contour plots at the 
same age. In the quadrant corresponding to radiation- 
confined shells, the H/3/H/3o line ratio is 0.35. This is 
consistent with our prediction, in Appendix [B] that /j on 
falls to ~ 0.3 in the limit \fr — > 0, il — > for the values of 
7 and /3 which characterize a 2 Myr cluster. 

6. DISCUSSION: IMPLICATIONS OF W bs SATURATION IN 
STARBURST GALAXIES 

On the basis of the theory developed above, we expect 
to find that individual H II regions (or regions within 
them) display line ratios consistent with ionization pa- 
rameters ranging up to the maximum value correspond- 
ing to radiation confinement. Although the same phys- 
ical process limits the inferred ionization parameter for 
observations on galactic scales, the value of the upper 
limit is not the same, because of the different measures 
used: usually geometric estimates (our U gcom ) in resolved 
regions, as opposed to pure line ratios (our U b s , or more 
precisely (U h s ) ga i) in unresolved ones. Moreover, aver- 
aging over many regions tends to blend together those at 
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^obs,max with those below it, and the selective dimming 
ol high-W bs regions by internal dust absorption favors 
those below the maximum. 

We estimate the upper limit of U gcom to be roughly 
10 _1 using the solutions of Appendix [B] by replacing 
^i,obs with the emission- weighted density n om . This is 
only approximate, as it does not account for the forma- 
tion of the lines used to derive n^bs- For unresolved re- 
gions we apply W bs,max ; in which we do account for line 
formation: because most of the observations employ the 
the [Ne III] /[Nell] line ratio of figure [51 we apply the up- 
per limit of 10 -1 ' 97 we found there. Part of the difference 
between the two upper limits can be seen in the internal 
structure of radiation-confined zones, in which S gcom is 
about 0.3 dex higher than the recombination-weighted 
average Si. om - The rest of the difference is presumably 
due to the details of line formation. 

To test these expectations against observed regions, we 
plot in FigureQ2]a sample of observationally-derived ion- 
ization parameters from the literature, divided according 
to the scope of the observations involved and according 
to the method used. As expected, the data range up 
to values which are approximately our estimated max- 
imum value for each method. For individually-resolved 
regions, the observed upper limit appears to be some- 
what higher than our prediction, but we attribute this 
to the approximate nature of our W gcomjmax (rather than 
to magnetic fields, say). For galactic-scale observations 
highest U hs is slightly lower than which 
could be a product of the averaging process. 

These results in hand, we turn to the importance of 
wind pressure in H II regions. 

For a resolved, spherical H II region with well- 
constrained observations of S, n e , and U b s , it is pos- 
sible to estimate the wind parameter Q by reference to 
a plot like our Figure |8j With only U b s or U gcoln , one 
obtains an upper limit on fi. By this logic, we can state 
quite firmly that the wind parameter f2 is often small in 
resolved H II regions: otherwise, our example values of 
Wobs in the upper panel of Figure [T2l would be much less 
than the theoretical maximum. 

The same argument applies, if somewhat more weakly, 
to entire galaxies: for them, log 10 U b s ranges up to about 
—2.3. This is far higher than what would be possible if 
they were adiabatic wind bubbles (Eq. \W\i unless they 
are exceptionally weak. This can be seen, for instance, 
in Figure [SJ where we use the values of (Mobs} ga \ and the 
product of total S and inferred n iiQ b s to place several star- 
burst galaxies within the parameter space of individual 
H II regions. Because this is an upper limit on the actual 
values of Sn em within each galaxy's H II regions, and be- 
cause (W bs) ga i cannot be higher than the typical values 
of U bs in those regions, we infer that the characteristic 
values of log 10 f2 cannot exceed about 0.25. 

While wind-dominated bubbles are likely to exist, we 
find no evidence for them either in individual regions 
or on galactic scales. We conclude from this that the 
leakage of wind energy past the IF, either by advection 
or radiation, is significant. 

7. CONCLUSIONS 

Our primary motivation in this study has been to as- 
sess whether the ionization parameter can be used to 



measure the importance of radiation pressure in the in- 
flation of H II regions in external galaxies, especially star- 
burst galaxies. We have been interested in how well the 
observational estimate U b s reflects the radiation-to-gas 
pressure ratio 2 in local conditions and its average over 
galactic scales. Because S oc U in each parcel of gas, 
and because U controls the ionization state and there- 
fore the emitted spectrum, it is possible, in principle, 
to derive the role of photon momentum in the kine- 
matic feedback from young star clusters into the nearby 
dense neutral gas. (Although supernovae tend to inject 
more energy, they are t hought to coup le very poorly to 
dense, star- forming gas; IMa tzncr 2002]-) Because emis- 
sion line ratios reflect local conditions, this informa- 
tion would be complementar y to the global analysis of 
lAndrews fc Thompson! (|2011D . To this end we explored 
the properties of strongly radiation-confined layers using 
the analytical theory of Drll (Appendix IB")) . considered 
a host of physical effects and scenarios which could affect 
Wobs in real regions (§[3] and §2]), and engaged in a suite 
of numerical models of dusty H II regions in quasi-static 
force balance with an interior pressurized wind bubble 
(§ El). Our models, though not meant to probe the full 
parameter space of H II regions, are the first to explore 
a wide range of both the radiation force parameter (r or 
and the wind pressure parameter (f2). 
Two effects complicate our plan to measure radiation 
pressure feedback using the apparent ionization parame- 
ter. 

First is the fact that strongly radiation-dominated re- 
gions compress their ionized gas into thin layers whose 
internal pressure gr adient balances the radiation force 
(|Binette et al.l ll997). The ionization parameter recorded 
in any line ratio reflects the ratio of radiation to gas 
pressure where the lines are produced, and in radiation- 
confined regions U \, s saturates at a maximum value. The 
resulting values of log 10 W b s ,max depend on the lines in 
question (§ [5J) and parameters such as the dustiness of 
the region and the ionizing spectrum, but are roughly 
~ —1.5, the value for which S ~ 1. For the particular 
neon line ratio plotted in Figure [SJ the upper limit is 
actually —1.97. Only significant levels of ram pressure 
or magnetic pressure, or serious departures from force 
balance, allow for Li b s to exceed this putative maximum 
( ^3.3ti3.4|) . and this point may be useful in diagnosing the 
physical states of unresolved regions. Conversely any sig- 
nificant pressure on its inner boundary, such as the pres- 
sure due to shocked stellar winds, reduces 2 and U Q bs 
(as well as the dust optical depth to) below their maxi- 
mal values (§ 13. ip . Our limited compilation of ionization 
parameters in individual regions (figure 1121 top panels) 
is consistent with the proposition that U hs varies up to 
the ma ximum imposed by ra diation confinement; for in- 
stance, [Snijders~eFal. (2007) find log 10 W o bs = —1-57 in 
individual regions in the Antennae Galaxies. 

Another complicating effect is the suppression of the 
apparent ionization parameter on galactic scales by se- 
lective dust absorption of ionizing photons in the regions 
where U is highest. Gas in the local Universe is suffi- 
ciently dusty (7 ~ 10) that the dust optical depth ex- 
ceeds the radiation-to-gas pressure ratio 3, both of which 
vary proportionally to U: all radiation-confined regions 
(and some others) are dust-limited. Such dust-limited 
regions convert a large portion of the stellar ultraviolet 
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radiation into far infrared continuum, limiting the effi- 
ciency /; on of line production. For observations which 
span a wide range of local conditions, this biases the in- 
ferred estimate of U (which we call (U hs) aa \ m * ms case ) 
toward regions of lower < 1 and U h s - 

Within the toy model for galactic H II region popula- 
tions we present in £|4.1[ it appears this selection effect 
suppresses (W bs) ga i by a small amount (<0.05 dex) rel- 
ative to an unbiased (i.e., starlight-weighted) average of 
^obs, primarily because the line efficiency /; on gradually 
declines by only 0.5 dex (to a minimum of about 0.35) 
as radiation confinement becomes strong. Although this 
result depends on the assumptions in our toy model, it 
strongly suggests that (U b s ) gal can be corrected for se- 
lective absorption and used to infer the physical char- 
acteristics of the galactic H II region population. Chief 
among these would be the ratio of radiation and gas- 
pressure forces averaged over a galaxy's ionized zones, a 
quantity we plot (for the toy model) in the top panel of 
figure [5| 

It has long been appreciated that the dust-bounded 
state contributes to the suppression of fine-structure 
cooling li nes relative to the far-infrared c o ntinuum in 
ULIRG s (fVoitl fl992t IBottorff et afl fl998t lAbel et al.l 
2009b), and that recombination luminosity under- 
represents the ionizing luminosity in this state; 

these effects provide observational checks on the degree 
of dust saturation one might infer from recombination 
line ratios. 

The galaxies collected in the bottom panel of figure 
Q21 display an upper limit of log 10 (W bs) ga i — —2.3, with 
some uncertainty. This is very close to the upper limit 
of -1.97 we found for U ohs derived from a [Ne III] / [Ne II] 
line ratio in young H II regions (3S|), especially consid- 
ering the bias induced by dust attenuation. We infer, 
therefore, that radiation pressure strongly modifies the 
ionized zones in these galaxies. 

A result of dust-induced bias is that the upper limit of 
(^obs) ga i will be weakly anti-correlated with metallicity. 
To test this, one must disentangle variations in U from 
variations in the incident spectrum due to changes in 
stellar metallicity, and possible influences of metallicity 
on the stellar initial mass function (and more subtly, on 
stellar rotation and multiplicity). 

Moreover, our analysis of the radiation-confined state 
(Appendix [B]) suggests a change in the properties of 
radiation-confined H II regions with metallicity above 
and below about \og 10 (Z/Z Q ) = —1.20 corresponding 
to 7 = 1. For higher values of 7, as we have already 
mentioned, dust absorption is significant in all radiation- 
confined zones. Further, for 7 > 1 radiation-confined 
zones are compact (the density scale height is much 
smaller than the radius) and the Lyman-a line pres- 
sure is negligible. None of these is true for 7 < 1: re- 
gions with significant radiation pressure extend inward 
toward the central source, are supported by line pressure 
as well as gas pressure, and are not significantly attenu- 
ated by dust. Moreover, the wind pressure will typically 
be lower on account of the lower opacity of stellar atmo- 
spheres. We leave a full investigation of this transition to 
future research, but point out that many works on photo- 

5 Estimated assuming T oc Z~ 0,2 , e; oc Z -0 - 04 , and oc Z. 



ionization in active galactic nuclei have been carried out 
in the limit 7 = 0. 

We firmly rule out the hypothesis that resolved H II 
regions are predominantly the ionized shells of energy- 
conserving wffid__bjjtbl^ (as envisioned for giant H II 
regions bv lOev fe~C larke 1997]) ■ The pressure of shocked 
winds in such models lowers U well below observed val- 
ues. Although this discrepanc y can be remedied b y dec- 
imating the wind luminosity ( Dop ita et al.l 12005) . it is 
more natural to attribute the reduction in wind pressure 
to leakage of hot gas beyond the ionization front, either 
via direct outflow or radiative losses. It is preferable to 
emp loy models which explicitly ac count for such leak- 
age (jHarper-Clark k Murray! 120091 KM09). Such mod- 
els form the basis for populations of H II regions in the 
toy model of £14. H and in a future companion paper (Ver- 
dolini et al. in prep.), in which line emission by star- 
forming galaxies will be reconsidered in greater detail 
than is possible here. 

It is possible to interpret entire galaxies in terms of 
models of individual H II regions, as we discussed in $6j 
This must be done with caution, because observations 
on galactic scales must necessarily average over a wide 
range of conditions, and are biassed by dust absorption. 
Acknowledging this, one can nevertheless derive an ef- 
fective upper limit on the wind pressurization Q from 
the fact that a larger value of this parameter would limit 
(Wobs)g a i below what is observed. 

For observational and theoretical studies of the ioniza- 
tion parameter, an important point is that the method 
used to infer U a h s affects the result, leading to system- 
atic offsets between methods. Traditionally line ratios 
are compared with suites of simulations in which IA \, S 
represents the interior of an ionized layer (often one of 
uniform density), a definition which corresponds to our 
Wgeom so long as the layer is thin. This definition en- 
counters a couple difficulties. First, the incident flux is 
attenuated by recombination and dust within the layer 
so that the characteristic value of U is lower than the 
incident value (e.g., equation [16]). Second, ionized gas 
can only be confined to a thin layer by an inner, pressur- 
ized region of hot gas, or (when 7 > 1) by strong radia- 
tion pressure; however, radiation confinement introduces 
a pressure and density gradient within the confined layer, 
and wind confinement limits the maximum values of S 
and U well below unity and 0.1, respectively. For these 
reasons, uniform-slab calculations are physically incon- 
sistent at high values of U. It is preferable to use suites 
of simulations involving spherical H II regions in quasi- 
static force balance, but a single sequence of these ex- 
cludes the role of a wind pressurization parameter (our 
Q). 

As an alternate standard, we advocate using, for U h s , 
the value of U within a uniform mixture of gas photons 
for which the emitted line ratios match those observed. 
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Figure 6. Normalized pressure (top) and electron density (bot- 
tom) profiles of H II regions, showing the structure of radiation 
pressure dominated (solid) and wind bubble (dashed) regimes, 
which were illuminated by the ionizing spectrum of a 2 Myr cluster 
of 10 42 ergs -1 and 10 39 ergs -1 luminosity, respectively. A drop in 
temperature beyond the ionization front leads to an increase in 
density which outpaces the declining ionization fraction, causing 
an uptick in n e in the outermost zones. 
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Figure 8. Logarithm of ionization parameter (solid blue con- 
tours) as a function of log 10 *I/ and log 10 f2, for illumination by 
a 2 Myr-old cluster. The log^Wobs values are evaluated from the 
line ratio [Nelll];^ 5 5^ m /[NeII],\i2.8i£tm- Dashed black contours 
denote the logarithm of Sn em . H II region of the central 500 pc 
of M82, NGC 3256, and NGC 253, are marked by open diamond, 
circle, and square, respectively. The maximum value corresponds 

to Wobs.max = —1.97. 
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Figure 7. One-zone [Ne III]/[Ne II] line ratio at various values 
of log 10 W, for hydrogen density log 10 = 2,3,4,5 (blue dashed 
line, red crossed line, cyan line with diamonds, magenta line with 
circles, respectively). Metallicity and ionizing spectrum are held 
fixed: Z e and the spectrum of a 2 Myr cluster. The [Ne III] /[Nell] 
ratio is sensitive to U and insensitive to nn, falling only 37% for 
densities which approach the critical density of the [Nell] line. 
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Figure 9. Same as in Figurc[8] but for ionization by a star cluster 
of age 4 Myr. The log 10 W D b s values are evaluated from the line 
ratio [NeIII]xi5.5S/im/[NeII]Ai2.81/im- Symbols that mark H II 
regions of M82, NGC 3256, and NGC 253, are the same as in 
Figure [8] 
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Figure 10. As in Figure [8] for ionization by a star cluster of 6 
Myr age. Values of log 1() W b s are evaluated from the line ratio 
[NeIII]xi5.55/Ltm/[NeII] A i2.8i/im- 
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Figure 11. Solid blue contours mark /i on , calculated as the ratio 
of the H/3 line to H/3o, its value in the absence of dust, as a function 
of log 10 \E" and log 10 O. Dashed black contours label log 10 (5n em ) 
as in previous figures. The ionizing spectrum corresponds to a star 
cluster age of 2 Myr. 
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Figure 12. Samples of individual H II regions and galactic-scale 
regions, plotted and divided into groups according to how W orjs 
is estimated. The thick horizontal line divides individual regions 
(upper portion) and galactic-scale regions (lower portion). We fur- 
ther distinguish the individual regions into two categories, based on 
whether U is estimated geometrically (W gC om = S/ {iTrR^ ¥ ri i ohs c)) 
or inferred from observed line ratios. Vertical dashed lines indicate 
theoretical maximum values W gC om,max (Appendix [Bj estimated 
from max(Sj igGom )) and W b s ,max (from [Ne III] /[Nell], figure^, 
as well as the value of W orjs at which dust absorption becomes signif- 
icant in the sense that /i on = 0.5 in individual young regions with 
our fiducial p arameters (Figure Q . Open circles - single H II re- 
gions in M82, McCradv & Graham (2007); Large open circle - H II 
region M82-A1, Smith et al. (2006); Open squares — re solved H II 
regions in the Antennae Galaxies, Sniidcrs ct al. ( 2007 |) ; C y an an d 
black left-pointing triangle - NGC 3256, Thornley et al. (2000); 
Cyan and black diamond - NGC 253, 1'1'hornlev et al.l (|2000| ); Cyan 
and black circle- the central 500 pc starburst region of M82, 
w here the bar indicates the poss ible range of log 10 U a h s modeled 
bv lForst cr Schreiber ct al. ( 2003) ; crosses — recali br ated logT n ^obs 
values of the galaxies in |Thornlcy ct al. (2000). Sniidc rs et alj 
(2007) report a value which corresponds to our W ge om, so we apply 
the approximate upper limit on this quantity. 
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APPENDIX 

A. FORCE BALANCE AND PRESSURE EQUILIBRIUM 
The momentum equation for fluid flow is 

d 

— pv + V • (pvv + PI) = f rad + f mag 

where I is the unit tensor and f ra( j and f mag are the magnetic and radiation force per unit volume, respectively. As 
elsewhere in the paper, P (with no subscript) is the gas pressure. Our analysis is based on a quasi-static idealization 
of an H II region in which v — > and Jr — > 0, so this simplifies to a statement of force balance: 

VP — frad ~t~ fmag- 

As discussed bv lMihalas &: Mihalasl (|1984D . f ra d = V ■ P ra d where P ra d is the radiation pressure tensor. In spherical 
symmetry, f ra( j = — J^P ra d — (3P ra d — w ra d)/r for radiation pressure P ra d and energy density it ra d- Splitting the radiation 
into radially free-streaming and isotropic components P ra d = P ra d,fs + P ra d,is with energy densities u ra d = P ra d,fs + 
3P ra d,isj we have f ra d = — ^P ra d,is ~ r ~ 2 gf r2 Pad,fs- Our free-streaming radiation is meant to represent attenuated 
radiation from the central source, whose luminosity at r is L(r) so that P r ad,fs( r ) = L{r) / (47rr 2 c), whereas the isotropic 
radiation field represents trapped line emission, especially Lyman a, as well as the diffuse Lyman continuum radiation 
emitted during recombinations to the ground state. 

Inserting this and f mag = — VP mag + (B • V)B/(4-7r) into the force balance equation and writing only the radial 
component, 

d d d & 

^-P = -^Pad.is - r- 2 — r 2 P rad ,f s - ^Pnag + (B • V)B r . 

or or or or 

Writing P tot = P + P ra d + P m ag and integrating across a range r ln — > r- m + Ar, 



47rr 3 c J r . 47r 

This shows that the change in Ptot is bounded across any range of r. Since L(r) is non-increasing, the magnitude of 
the first integral is maximized if radiation and matter do not interact, so that L(r) remains constant: 



L(r) , , J rf„ \ A 

V / J„ ^ o D _ . f„ \ I 1 ml , A 



47rr 3 c 



dr < 2P rad ,f s (r in ) 1 - -f±- — > 4 — p. a d,f s (r in ) 



\r 



where the arrow represents the limit Ar <C r- m . Moreover, the magnetic tension term (B • V)P r /(47r) is of order 
2P m ag/Pc if Rc is the curvature radius of the magnetic field. When ionized gas is compressed in a thin shell of width 
Ar, R c likely to be intermediate between Ar and r. We can therefore say in general that when the shell is thin, 

|AP tot | < ArP rad .fs^in , ^ -fmag 



4— lau - la 111 +2- 



Ptot T Ptot Rc Pot 

which is typically quite small: force balance implies near-constancy of the total pressure. 

B. DUSTY, RADIATION-DOMINATED H II REGIONS: THE PLANAR LIMIT 

We are interested in the properties of the Drll solutions in the limit that ionized gas is restricted to a geometrically 
thin shell, as occurs in radiation-dominated regions when L > 1 and (as we see below) 7 > 1. We restrict equations 
(l)-(3) of Drll to the limit of a thin shell (constant r) and combine them to find 

% = -4>-~n (Bl) 
(we use r for the local dust optical depth, reserving 7^ for that of the entire region) and 

d:f) 



where n = nj 7i cri for a characteristic density 



r-= 1 - 1 ( ( f > + Pe- T +n) (B2) 



n ch = , ( B3 ) 



which is defined so that, when n — <f> — 1, dust absorption and ionization consume ionizing photons at equal rates. 
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The solution with n = and <fi — 1 at r<j = is 

*( T ) = 7(/3 + l)e-^-^-/3[(7-l)e- + l]-l (B4) 

7-1 

and 

n( T ) = I±£ [l - e -r(7-i)/7j (B5) 
7-1 L J 

valid up to t — Td, for which ^(r^) = 0, so that 

70S + l)e^ fcTd = /3( 7 - l)e- Td + /3 + 1 

for k — 7/(7— 1), which shows that when 7> 1 and k —> 1, Td —> ln[(7 + /?)/(/? + 1)]Q Because of the planar geometry 
(Appendix [XJ these solutions obey pressure equilibrium, 47rr 2 p = (L n e~ T + Li<p)/c, or in dimensionless terms, 

7 n + + /3e~ T = 1 + j8 (B6) 

so that, at the IF, ra = [1 + (1 - e~ rd )/3]/7. 
In the special case 7=1, we find 



Within these solutions, 



(r) = (l + /3)(l-r)-/3e- T ; h(r) = (1 + 0)r. (B7) 



7 \l d ndT. 



(B8) 



and the fraction of ionizing starlight caught by gas is 



/ion = 1 - /dust = / ndr 







1 + /3 



(B9) 



7 ^ 

r.i — — 

7-1 

~1 - cxp[0.7(lnr d - 1) - (lnr d ) 2 /6], (BIO) 

where the final approximation is correct to 2% for the entire range of figure IT3l In the very dusty limit 7 » 1, for 
which Td ~ 3, the recombination luminosity is suppressed by a constant fraction which is roughly (1 + /3)T<j/(7 — 1). 

The radial coordinate can be determined by integrating dr = dr/(nad). Because of the low density at the inside of 
the slab, radiation force is transmitted to the gas only through the dust, leading to a constant acceleration and (since 
the temperature is assumed to be c onstant) an exponential interior density distribution n cx exp[(l + P)nch^d r / 1\- 

Note, however, that equation |B5] implies a negative density if 7 < 1; this simply means that the solution is incompat- 
ible with the imposed inner boundary of zero density and optical depth. For such low dust opacity, the solution must 
extend inward to join the spherical solutions of Drll, or meet an inner, pressurized region of hot gas. The difference 
between compact (7 > 1) and extended (7 < 1) radiation-confined regions (with no inner hot gas) is apparent in the 
top panels of Drll's figure 2. For 7 < 1, radiation-confined regions with r > 1 are still concentrated toward the 
ionization front, but the scaling of radiative acceleration with the neutral fraction allows gas to permeate the interior. 
This effec t would only be ac centuated by the pressure of Lyman a line photons, which becomes significant when dust 
is scarce (| Arthur et alJl2004f >. 

The net dust optical depth Td and ionization efficiency /i on , which reach their maximum and minimum values 
(respectively) in the radiation-confined state, are plotted in figure [TBI 

In figure [T4l we plot two characteristic values of the ionizing radiation-to-gas ratio in the radiation-confined state: 
the recombination- weighted average Si jCm , and the geometrical Si igeom . The latter is calculated using the emission- 
weighted density n i cm . Estimates of Z// ODS and U geom derived from these limits are only approximate, because line 
emission is not strictly proportional to the recombination rate. 
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Figure 13. Dust optical depth, emission-averaged ionizing radiatio n prcssur e-to- gas ratio, and fraction of ionizing radiation absorbed by 
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recombination lines, constant temperature, identical dust cross sections for non-ionizing and ionizing radiation (whose luminosity ratio is 
/3), and no pressurized inner region (f2 = 0). The parameter 7 indicates the relative importance of dust in the consumption of ionizing 
photons in a state of equal gas and radiation pressures. 
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However the latter is better estimated using the models of [[5] for the specific lines in question. 
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